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Abstract 

In this paper we present an analytical study of a subgrid scale turbulence model of the three- 
dimensional magnetohydrodynamic (MHD) equations, inspired by the Navier-Stokes-a (also known 
as the viscous Camassa-Holm equations or the Lagrangian-averaged Navier-Stokes-a model. Specifi- 
cally, we show the global well-posedness and regularity of solutions of a certain MHD-a model (which 
is a particular case of the Lagrangian averaged magnetohydrodynamic-a model without enhancing 
the dissipation for the magnetic field). We also introduce other subgrid scale turbulence models, 
inspired by the Leray-a and the modified Leray-o? models of turbulence. Finally, we discuss the 
relation of the MHD-a model to the MHD equations by proving a convergence theorem, that is, as 
the length scale a tends to zero, a subsequence of solutions of the MHD-a equations converges to a 
certain solution (a Leray-Hopf solution) of the three-dimensional MHD equations. 
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1 Introduction 

We consider the three-dimensional magnetohydrodynamic (MHD) equations for a homogeneous incom- 
pressible resistive viscous fluid subjected to a Lorentz force due to the presence of a magnetic field. 
The MHD involves coupling Maxwell's equations governing the magnetic field and the Navier-Stokes 
equations (NSE) governing the fluid motion. The system has the form 

dv f 

-^ + {v-V)v- i^Av + Vtt + - V|Bp = (B • V) B, 

^ + (v-V)B-(B-V)v-'nAB = 0, ^^-^^ 
ot 

V - v ■ B = {), 

where v (x, t), the fluid velocity field, B (x, t), the magnetic field and tt, the pressure, are the unknowns; 
> is the constant kinematic viscosity and > is the constant magnetic diffusivity. 
Current scientific methods and tools are unable to compute the turbulent behavior of three-dimensional 
(3D) fluids and magnetofluids analytically or via direct numerical simulation due to the large range of 
scales of motion that need to be resolved when the Reynolds number is high. For many purposes, it might 
be adequate to compute only certain statistical features of the physical phenomenon of turbulence and 
much effort is being made to produce reliable turbulence models that parameterize the average effects of 
the fluctuations on the averages, without calculating the former explicitly. Motivated by the remarkable 
performance of the Navier-Stokes-a (NS-a) (also known as the viscous Camassa-Holm equations (VCHE) 
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or the Lagrangian-avcragcd Navicr-Stokcs-Q (LANS-a)) as a closure; model of turbidcncc in infinite chan- 
nels and pipes, whose solutions give excellent agreement with empirical data for a wide range of large 
Reynolds numbers [3-5], the alpha subgrid scale models of turbulence have been extensively studied in 
recent years (see, e.g., [3-7,14,15,24,25,30-32,39,48]). 

A justification of the inviscid NS-a model can be found, for example, in [5,8,21,22,35]. 

An extension of the NS-a model to the nondissipative MHD is given, e.g., in [20]. The model 
was obtained from variational principles by modifying the Hamiltonian associated with the ideal MHD 
equations subject to the incompressibility constraint. Then the dissipation is introduced in an ad hoc 
fashion in analogy to the NS-a, following [3-5, 15]. Specifically, the flow Lagrangian of the ideal MHD is 
given by 



/:[u,D,B] = I (^^D\u\'-n{D-l)-^\B\^^ 



dx 



with volume preservation for the pressure. Here the volume element D{x,t) = (det (dX/da) {a,t))~ at 

X = X{a,t), where X{a,t) is the Lagrangian fluid trajectory, dX /dt{a,t) = u{x,t) (sec [18]). First, the 
Lagrangian is averaged and approximated using a form of Taylors hypothesis (see, e.g., [19]) to obtain 



dx. 



then the Hamiltonian principle is applied (see, e.g., [22]) to produce an ideal MHD-a model (eq. (1.2) 

with h' — 7] = 0). Adding viscosity and diffusivity provides the MHD-a (or the Lagrangian-averaged 
magnetohydrodynamic-a (LAMHD-a)) model 

dv ^ 3 

-^ + {u-W)v + ^ VjVuj - uAv + Vp + ^{Bs)jVBj = {Bs ■ V) B, 

^ + {u-y)Bs-{Bs-y)u-r]AB = 0, (1-2) 

v={l- a^A) u, -B = (1 - aifA) Bs, 
Vtx = V- v = VB« = VS = 0, 

where u and Bs represent the unknown 'filtered' fiuid velocity and magnetic fields, respectively, p is the 
unknown 'filtered' pressure, and a > 0, aM > are lengthscalc parameters that represent the width 
of the filters. At the limit a = 0, aM — 0, we formally obtain the three-dimensional MHD equations. 
The LAMHD-a model was investigated numerically in periodic boundary conditions in two [38,41] and 
three [37] space dimensions against direct numerical simulations. In [41] the Karman-Howarth theorem 
was extended to LAMHD-a equations. The LAMHD-a model was also studied in [26] in the context of 
convection-driven plane layer geodynamo models. 

We tend to think about the a models as a numerical regularization of the underlying equation, which 
makes the nonlinearity nulder, and hence the solutions of the modified equation are smoother. This 
is contrary to the hyperviscosity regularization [33] and nonlinear viscosity [28,29,43], which lead to 
unnecessary extra dissipation of the energy of the system. To emphasize this numerical analysis point of 
view, we observe that recently a Leray-a model of the inviscid Burgers equation 

dv dv „ ,^ 

which is 

9t ^ " dx n 4) 

has been introduced in [1] and [45]. Regular unique solutions of (1.4) exist globally and it was shown 
computationally in [1] that the solutions of (1.4) converge to the unique entropy weak solution (see, 
e.g., [40,44,45]) of (1.3). Notice that there is no dissipation in (1.4), and the norm of v" is preserved. 
On the other hand, the viscous regularizing approach, which is usually taken for the Burgers equation, is 
achieved by introducing an artificial viscosity term in (1.3) and obtaining the viscous Burgers equation 
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This model gives a smooth sohition , which c;onverges in the appropriate norms to the unique entropy 
weak solution (see, e.g., [40]). However, the energy of is decaying in time at a much higher rate than 
the decay expected for the entropy weak solution. Hence, the advantage of introducing the Leray-a 
model (1.4) for Burgers equation. This simple example clarifies our numerical approach of why we insist 
on making the nonlinearity milder instead of adding additional viscous or hyperviscous terms. This 
approach has been discussed further in [2] in the context of Euler and Navier-Stokes equations. 

Filtering the magnetic field, as it is done in [26,37,38], is equivalent to introducing hyperdiffusivity 
for the filtered magnetic field Bg, due to the term — rya^A^Bg in (1.2), which we think is unnecessary. 
Taking the numerical analysis point of view discussed above we prove the well-posedness of a certain 
MHD-a model without introducing extra dissipation for the magnetic field, i.e. we filter only the velocity 
field, but not the magnetic field and obtain the following regularizing system of (1.1) 

uAv + Vp + iv|B|2 = (S • V) B, 

r]AB = 0, (1-6) 



instead of the system (1.2). 

As a models are some sort of regularizing numerical schemes, we would like to make sure that they 
inherit some of the original properties of the 3D MHD equations. Formally, three ideal, i.e. v = t] = 0, 
quadratic invariants of the system (1.6) could be identified with the invariants of the original ideal 
3D MHD equations under suitable boundary conditions, for instance, in rectangular periodic boundary 
conditions or in the whole space R^. Namely, the energy E" = ^ Jq{v (x) ■ u{x) + \B{x)\'^) dx, the cross 
helicity iJg — ^ In '"i^) ' B{x)dx, and the magnetic helicity H'^ — k In ■^i-"') ' B{x)dx, where A is the 
vector potential, so that B = 'V x A; and they reduce, as a ^ 0, to the ideal invariants of the MHD 
equations. 

There are other possible alpha subgrid scale models that can be shown to have global existence and 
uniqueness. For instance, inspired by the Leray-a [6, 7, 17, 23, 48] and modified Leray-a (ML-a) [25] 
models of turbulence, we formulate similar MHD alpha models, we refer to them as Leray-a-MHD and 
ML-a-MHD models, respectively. The Leray-a and ML-a models of turbulence reduce to the same 
closure model for the Reynolds averaged Navier-Stokes equations in turbulent channels and pipes as 
the NS-a model under the corresponding symmetries [6,7,25], which, as we mentioned above, compares 
successfully with experimental data for a wide range of Reynolds numbers. This comparison means that 
the Leray-a and the ML-a models as well as NS-a equations could be equally used as subgrid scale 
models of turbulence. 

Specifically, we consider the following version of the three-dimensional Leray-a-MHD model 

-^ + {u-V)v- uAv + yp+ - V|B|2 = (B . V) B, 
OB 

— + (w • V) B - (B ■ V) v - ?7AB = 0, (1.7) 

v= {1- a^A) u, 

V-ix = V- i; = V- B = 0. 

Formally, the term (B • V) u comes from requiring in the ideal (v = = r]) case the conservation of 
energy i?" ~ ^ J^^ (\v{x)\'^ + |B(a;)|^) dx (under suitable boundary conditions). While the requirement 
for the system to have an ideal invariant corresponding to the cross helicity Hq = ^ IqV{x) ■ B{x)dx 
leads to the term (u- V)B. Contrary to the MHD-a model (1.6), where we establish the existence 
and uniqueness, for the 3D Leray-a-MHD model (1.7) we are able to establish only the existence of 
weak solutions, as in the case for the original MHD equations (1.1). In this case, the term (B • V) v is 
problematic as in the usual 3D NSE and MHD. However, in the two dimensional case the existence and 
uniqueness of weak solutions can be shown (similarly to the proof given for the model (1.6) in section 3) 



dv 
'dt ■ 

dB 

'dt 

v= {1- Q^A) u, a > 0, 
V-u = V-t; = V-B = 0, 



\- {u ■V)v + Vj Vuj — 
+ {u-V)B-{B-V)u- 
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for the following 2D-Leray-Q!-MHD model 



— + {u-V}v- i^Av + Vp + - V|Bp = (S • V) B, 
QB 

— + V)B-(B-V)u-r/AB = 0, (1.8) 

v= {1- a^A) u, 

V ■u = V -v = V ■ B = 0. 

For this system, due to the identity (it • Vu) • Ait — (for the periodic 2D case and divergence free u), 
the ideal invariant corresponding to the energy is = ^ Jq{v {x) u{x) + |B(x)p) dx. At the moment 
we are unable to find a conserved quantity in the ideal version of (1.8) that can be identified with a 
cross helicity. The mean-square magnetic potential, given by ^4. = 5 /q \ip{x)\'^dx, where B = V^Vj is 
conserved in the ideal case. We note that it appears that there is no conserve quantity that could be 
identified with energy for the 3D version of (1.8). 

The three-dimensional Modified-Leray-a-MHD model, for which the well-posedness can be proved in 
a similar way as for the model (1.6), is given by 

^ + {v-V)u- uAv + Vp+ ^ V|B|2 = {B-V) B, 
dB 

— + iu-V)B-iB-V)u-7]AB = 0, (1.9) 

V = {1- a^A) u, 

\/ .u = V ■v = V ■ B = 0, 

where the term (B • V) tx comes from requiring the conservation of energy (in the ideal case, with 
periodic boundary conditions or in M^) E" = ^ Jq{v (x) u{x) + |B(x)p) dx. Also, the system conserves 
the magnetic helicity H'^ = i J^-^A{x) ■ B(x)dx. At the moment we are unable to find a conserved 
quantity in the ideal version of (1.9) which can be identified with a cross helicity. 

The main goal of this paper is to establish the global existence, uniqueness and regularity of solutions 
of the three-dimensional MHD-O! equations (1.6) subject to periodic boundary conditions (similar results 
also hold in R^). We emphasize again that we consider a version of the MHD alpha models, where only 
the velocity field is filtered, while the magnetic field remains unfiltered. We note that in the case of 
filtcaing the magnetic field, as in (1.2), one has hyperdiffusivity for the filtered magnetic field Bg and 
the proof of the existence and uniqueness of regular solutions of (1.2) is deduced in a similar way. 

We start by introducing some preliminary background and the functional setting in section 2. In 
section 3 we show the global well-posedness of the MHD-q: subgrid scale model of turbulence (1.6). We 
remark that using the Gevrey regularity techniques developed in [16] (see also [13]) one can show that 
the solution of the MHD-a model becomes instantaneously analytic in space and time. As a result of 
this Gevrey regularity, one deduces the existence of a dissipation range in the energy spectrum in which 
the energy decays exponentially fast as a function of the wavenumber, for k larger than the dissipation 
length scale (see [11]). One can also establish, in the forced case, the existence of a finite dimensional 
global attractor, a subject of future work. In section 4 we relate the solutions of the MHD-O! equations 
to those of the 3D MHD as the length scale a tends to zero. Specifically, we prove that one can extract 
subsequences of weak solutions of the MHD-a equations which converge as a ^ 0+ (in the appropriate 
sense) to a Leray-Hopf weak solution of the three-dimensional MHD equations (1.1) on any time interval 
[0, T], which satisfies the energy inequality 

\v it) |2 + \B (t) [^ + 2 /* (Hl^(s)f + v\\B{s)f) ds < \v (to) f + \B (to) f 

for almost every to £ [0,T] and all t G [^Oj?^]- Also, if the initial data is smooth a subsequence of 
solutions converges for a short interval of time, that depends on the initial data, v, rj and the domain, 
to the unique strong solution of the MHD equations on this interval. Thus the a models can be viewed 
as a regularizing numerical method. Section 5 contains a discussion summarizing our results. 



4 



2 Functional Setting and Preliminaries 

Let ft be the L-periodic three-dimensional box fl = [0,L]^. We consider the following MHD-a subgrid 
scale turbulence model, which we introduced in (1.6), subject to periodic boundary condition with a 
basic domain f2, 



^ +{u-V)v + J2 Vttj- - ^^■f' + + ^V|-Bp = (-B • V) B, (2.1a) 



dv 
'dt 

— + {u-V)B -{B ■S/)u-r]AB = 0, (2.1b) 

v={l- a^A) u, (2.1c) 

V-w = V-i; = V-B = 0, (2. Id) 

u{x,0) =u"\x), (2.1c) 

B{x,0) = B'"{x), (2.1f) 

where u represents the unknown 'filtered' fluid velocity vector, p is the imknown 'filtered' pressure, and 
a > is a lengthscale parameter which represents the width of the filter. At the limit a = we formally 
obtain the three-dimensional MHD equations (1.1), where u is the Eulerian velocity field and p — ^\u\^ 
is the pressure. Notice that we chose to smooth only the velocity field and not the magnetic field, thus 
we do not introduce hyperdiffusivity for the magnetic field, as it is for the filtered magnetic field in (1.2). 
We consider initial values with zero spatial means, i.e., we assume that 

/ ti"da;= [ B'"dx = 0, (2.2) 
Jn Jn 

then from (2.1a) and (2.1b), after integration by parts, using the spatial periodicity of the solution and 
the divergence free condition (2. Id) we have {d/dt) vdx = 0, {d/dt) Bdx = and (d/dt) udx = 0. 
Namely, the spatial mean of the solution is invariant under time. Hence, by (2.2), ^^vdx = J^udx = 
LBdx = Q. 

Next, we introduce some notation and background following the mathematical theory of NSEs, see, 
for instance, [9,29,46,47]. Let Lp{^) and H"^{n) denote the Lebesgue spaces and Sobolev spaces 
respectively. We denote by | • | the L^-norm, and by (■, ■) the L^-inner product. Let X be a linear subspace 
of integrable functions defined on the domain i7, we define X := {ip ^ X : (p{x)dx = 0} and V = {(p : 
if is a vector valued trigonometric polynomial defined on fl, such that V • <^ = and (p{x)dx = 0}. 
The spaces H and V arc the closures of V in L^{il) and in H^{n) respectively; observe that H^, the 
orthogonal complement of H in L^(f7) is {Vp : p E H^{^1)}. Let P<j : (fi) H he the Hclmholtz-Leray 
projection, and A = —Po-A be the Stokes operator with domain D{A) = {H'^{^1) DV). In the periodic 
boundary conditions A = — A|£i(^) is a self-adjoint positive operator with compact inverse. Hence the 
space H has an orthonormal basis {wjyj^Li of cigenfunctions of A, Awj = \jWj, with < Ai < A2 < . . ., 
Xj ~ j2/d^-2^ ggg^ g_g^^ [9^3g]_ One can show that V = D {A^/^). We denote ((•, •)) = (A^/^-, A^a.) and 
II • II = |A^/^ • I the inner product and the norm on V, respectively. 

Following the notation of the Navier-Stokes equations and those of [15], we denote 

B {u, v) = Pa [{u • V) w] , U,V GV, 
B {u, v) = Pa [(V X v) X u] , u,v & V. 

Notice that 

[B [u, v) ,w) = — {B {u, w) ,v) , u,v,w G V, 

and due to the identity 

(6 • V) a + ^ ajVbj = -6 x (V x a) + V (a ■ 6) , (2.3) 



^B {u, v) , = {B (it, v) ,w) — {B (w, v) ,u) . 
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The definitions of B {u, v) and B {u. v) and tlic above algebraic identities may be extended to larger 
spaces by the density of V in the appropriate space each time the corresponding trilinear forms are 
continuous. The extensions of the bilinear forms B and B (which we also denote B and B) have the 
following properties 

Lemma 2.1. 

(i) Let X he either B or B. The operator X can be extended continuously from V xV with values in 
V (the dual space ofV). In particular, for every u,v,w £ V, 



Moreover, 

which in turn implies that 
Also 

and hence 



\{X{u,v),w)y,\<c\unu\\^/^v\\\\w\\. 

{B {u, v) ,w) = — {B {u, w) ,v) , u,v,w G V, 
{B {u,v) ,v) = 0, u,veV. 
(^B {u, v) , = {B {u, v) ,w) — {B {w, v) ,u) , u,v,w G V, 
^B {u, v),uj =0, u,v e V. 



(a) Furthermore, let u G D{A), v € V,w G H and let X be either B or B then 

\{X{u,v),w)\<c\\u\\'/^\Aunv\\\w\. 

(Hi) Letu gV,v G D{A), w G H then 

\{B {u,v) ,w)\ <c||w||||t; \W\. 

(iv) Let u G D{A),v G H, w gV, then 

\{B{u,v),w)v'\ < c||M||i/2|Ati|i/2 

(v) Letu,v,w G V, then 



< c||m||||v|||w|^/^||w||^/^. 



{B {u,v),w)v 

(vi) Let ugH,vgV, wgD {A} and let X be either B or B then 



{X {u,v),w)jj^^y 
(vii) Let u G V, V G H, w G D (A) then 



< c\u\\\v\\\\w\\^/^\Aw\'/^. 



B (u,v) ,w) < c 
/D{Ay 

and hence by Poincare inequality. 



B {u,v) ,w 



(via) Let u G D {A) , v G H, w G V then 



DiA)' 



(2.4) 

(2.5) 
(2.6) 
(2.7) 
(2.8) 

(2.9) 

(2.10) 

(2.11) 
(2.12) 

(2.13) 



(2.15) 



B{U,V),W^^^\ < C (lluf /2|Au|V2|„|||^|| + |yl„||t,||t«|l/2||^t,||l/2^ . (2.16) 



In this lemma and throughout the paper c denotes a generic scale invariant constant. 
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Proof. The proof of (i) can be found, for example, in [9,46,47] for B and in [15, Lemma l(iii)] for B. 
To prove (ii) we first consider the case where u,v,w £V 



\{B {u,v),w)\ = 
\{B {u,v),w)\ = 



• wdx 



/ {u-V)v 

Jn 

/ u X {V X v) 
Jn 



wdx 



hence 



\{X{u,v),w)\<c\\u\\^^ \\Vv\\^, \\w\\^,. 
By Agmon's inequahty in three-dimensional space, see, e.g., [9], 



we obtain 



ll</'llLoo<||</'||^^f||</'||^4' 



\iX iu,v),w)\ < c\\u\\^^^\Au\^/^\\v\\\w\. 



Since V is dense in D{A), V and H we conclude the proof of (ii). 

To prove (iii) we recall the following Sobolev-Ladyzhenskaya inequalities (see, e.g., [9,29]) in 3D 



lL6 



<C 



|HL3<C|<^|1/2||</,||V2, 



for (A € V. Then we have 



{B {u,v),w) 



In 
< c\\u 



/ (it • V)v ■ wdx 

Ja 



w 



L2 



<c||u|||Vt;|i/2||Vt;||i/2|u;| 
<cM\\vr/'\Av\'/'\w\. 

The proof of (iv) is a direct result of the (ii) due to the symmetry (2.5). The proof of (v), (vi), (vii), 
(viii) can be found in [15, Lemma 1 (iii,w,v,vi)]. □ 

Using the above notations and the identity (2.3) we apply to (2.1) to obtain, as for the case of 
the NSE, the equivalent system of equations (see, e.g., [46] and [12]) 



— + B{u,v)+ vAv = B {B, B) , 

^ + B{u,B)-B {B, u) + r]AB = 0, 

u{0) = u*", 
B(0) = B*". 



(2.17a) 

(2.17b) 

(2.17c) 
(2.17d) 



Definition 2.2. Let T > 0. A weak solution of (2.17) in the interval [0,T], given u{0) = u'" G F (or 
equivalently v'" G V) and B (0) = B*" e iJ , is a pair of functions u, B, such that 

ueC ([0, T];V)n ([0, T] ; D (A)) with ^ e ([0, T] ; H) 

(or equivalently veC ([0, T] ; V) D L"^ ([0, T] ; H) with ^ G ([0, T] ; D {A)')) and 

B G C ([0, T] ; if ) n ([0, T] ; V) with ^ G ([0, T] ; F') , 
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satisfying 

^v,w) +(b{u,v),w) +u{v,Aw) = {B{B,B),w)y,, (2.18a) 

/ D{AY ^ ' ^'-^> 

j^B, 1^^^ + {B {u, B),0- {B {B, u),^)+v {{B, ^)) = (2.18b) 

for every w G D (A) , ^ gV and for almost every t G [0, T]. 

The equation (2.18) is understood in the following sense: for almost every to,t € [0,T] 

{v{t),w)-{v{to),w)+ [ (b {u{s) ,v{s)) ,w) ds + iy [ {v{s),Aw)ds (2.19a) 

Jto ^ /-D(^)' Jto 

= [ {B{B{s),B{s)),w)y,ds, 

Jto 



(B (t) ,^)-iB (to) ,^)+ {B{u {s) , B (s)) ,Ods- {B [B {s) , u (s)) , $) ds (2.19b) 

Jto Jto 



+ V [\{B{s),^))ds = 0. 

Jto 



When ti™ e D{A) (or equivalently € H) and G F we call a strong solution of (2.17) in the 
interval [0, T] the solution that satisfies 

BeC ([0, T];V)n ([0, T] ; D(A)) , « G C ([0, T] ; D(A)) n ^^([o, T] ; DiA^/^)) 

(or equivalently t; G C ([0, T];H)n ([0, T] ; D{V))). 



3 Global existence and uniqueness 

In this section we show the global well-posedness of the MHD-a model (2.1) or equivalently (2.17). 

Theorem 3.1. Let w'" G V, B'" G H. Then for any T > there exists a unique weak solution u, B of 
(2.17) on [0,T]. Moreover, this solution satisfies 

as well as the energy equality 

\u (t) \' + a'\\u (t) f + \B (t) \' + 2 f (Kll«(s)f + a'lAu^sr) + r,||B(s)f ) ds 

Jto 

= |«(io)|' + a'||M(io)ll' + l^(io)l', 0<to<t<T. (3.1) 

We use the Galerkin approximation scheme to prove the global existence and to establish the necessary 
a priori estimates. Let be an orthonormal basis of H consisting of eigenfunctions of the operator 

A. Denote = spanjwi, . . . , Wm} and let be the i^-orthogonal projection from H onto iJ^- The 
Galerkin approximation of (2.17) is the ordinary differential system 

+ PmB {Um, Vm) + yAVm = PmB {Bm, Bm) (3.2a) 

dB 

+ PmB {Um, Bm) - PmB {Bm, Um) + r]ABm = (3.2b) 

Vm = Um+ a^AUm (3.2c) 

Um (0) = P^m'" (3.2d) 
Bm (0) = PmB'^. (3.2e) 

Since the nonlinear terms are quadratic, hence locally Lipschitz, then by the classical theory of ordinary 

differential equations, system (3.2) has a imiquc solution for a short interval of time (— r^, Tm)- Our goal 
is to show that the solutions of (3.2) remains finite for all positive times, which implies that Tm = oo. 
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3.1 ff^-Estimate of Um, L^-Estimate of 

We take the inner product of (3.2a) with Um and the inner product of (3.2b) with and use 
(2.6), (2.8), (2.5) to obtain 

~ {\Umf + a'^WUmf) + V (||Mmf + c?\AuJ^) = [B (B„, B„) , Um) , (3.3a) 
~^\Bm? + ri\\B^f = -{B{Bm,Bm),Um). (3.3b) 
Now, by summing up (3.3a) and (3.3b), we have 

(|«™|^ + a"||«™f + \Bm?) + V {\\um\? + o?\Aum?) + vWBmf = 0. (3.4) 
We denote ji = min {v, rf\ and obtain 

\jt ^ "'ll'"™"^ + 1^™!') + dl"*™!!' + (^Mum? + \\Bmf) < 0. (3.5) 

Using Poicare's inequaUty we get 

^ (|w„|2 + a^Umf + |-B„|2) + 2/iAi (|u„|2 + a2||u^||2 + |_B„|2) < 0. 

and then by Gronwall's inequaUty we obtain 

\um it) P + a'Wum it) ||2 + \Bm it) P < e-^/^^i* (|m„ (0) ^ + a^jj^^ (q) f + |b„ (0) p) . 

Hence 

\um it) P + a^Wum it) f + |B„ it) |2 < fci := + a^ll^-f + (3.6) 

for all t > 0. 

This implies that Tm = oo. Indeed, consider [0,T™°^), the maximal interval of existence. Either 
jroax _ g^jj^^ g^j.g j^Qng^ or T™"*^ < (X) and we have hmsup(^(2.„ax)- {\um it) ^ + \Bm it) P) = oo, 
a contradiction to (3.6). Hence we have global existence of Um, Bm, and hereafter we take an arbitrary 
interval [0,T]. 

Integrating (3.4) over the interval (s, t) and using the estimate (3.6) we obtain that, for all < s < 

2 f {vi\\Um{T)f+a''\AUm{T)\^)+r)\\BmiT)\\^)dT<k^. (3.7) 
J s 

3.2 ff^-Estimate of u^, ff^-Estimate of 

By taking the inner product of (3.2a) with Au^ and the inner product of (3.2b) with AB^ we have 

ld_ 
2dt 

1 d 



\Umf + a^\AUm\^) + ly (\AUm\'^ + a^\A^/^Um\^^ = iBiBm,Bm),AUm) - (^BiUm,Vm),AUr, 



(3.8a) 



~\\Bm\\'' + v\ABm\'' = {BiBm,Um),ABm)-iBium,Bm),ABm). (3.8b) 
First, we estimate the nonlinear terms. By (2.12) we have 

(^BiUm,Vm),AUm)\<c{X^' +a'') \\Um\\\AUm\'^^\A^/^Umf/\ (3.9) 

To bound the term |(_B (£,„, B„i) , Aum)\ wc use (2.11) 

\iBiBm,Bm),Aum)\<c\\Bm\\'/^\ABm\'/^\Bm\\A^/^Um\. (3.10) 

By (2.9) we have 

\{B iBm,Um) ,ABm)\ < c\\Bm\\'/^\\um\\\ABmf/^ (3.11) 



and by (2.10) 

\{B{Ura,B^),AB^m)\<c\\Bra\\^/^\\Um\\\ABmf/\ (3.12) 

Now, summing up (3.8a) and (3.8b), we obtain 

= {B {Bm,Bm) ,AUm) - (-B {Um,Vm) ,AUmj + {B (-B^ , Um) , AB„) - {B {Um,Bm) ,ABm) ■ 

(3.13) 

By (3.9), (3.10), (3.11) and (3.12) and several applications of Young's inequality we reach 

^ (llw^f + a'\Au^\^ + IIB^f ) + {\Aum\^ + a^lA^/^^^l") + r;|AB„|" 

< c(a2i/)-3 (X-^ + a2)4 ||„„||4|^„^j2 ^ c(a2^)-2,y-i||B„||2|B„|4 + cr,''\\Bmf\\um\\\ (3.14) 
Integrating over (s,t) and using (3.6), (3.7) we obtain 

\\um it) ir + a^lAurn it) |2 + [t) f + (z. (|Aw,„ (r) p + a2|A3/2M„ (r) l') + 771 AB„ (r) |') dr 

< \\um is) f + a^\Aum is) ^ + (s) f + Ku (3.15) 

where we denote 

(i) Now, if u'" G r> (A), B'" e V, we have 

(Of +a'l^Mm (t)P + ||B™ (t)f 

+ y (i. (|Atx„ (r) p + a^l^^/^tx™ (r) |') + v\AB^ (r) |') dr 

< ||ti'"||2 + a2|ylu'"|2 + ||B'"f + :=fc2. (3.16) 

(ii) Otherwise, if m*" ^ (^), S'" ^ wc integrate 

it) f + a^\Aum it) |2 + \\Bm it) f < \\um (s) f + a^\Aum is) P + is) f + i^l 
with respect to s over (0, t) and use (3.7) to obtain 

t (||m„ it) ||2 + a^\Au„, it) p + ||B„ it) f ) < i-fci + Kit, 



hence for f > 



||tx„. (t) f + a^\Aum it) P + \\B^ it) f <K, + ^/^"'fci := ^2 (i), (3.17) 



and thus 



J (u (\Aum (r) |2 + a2|A3/2w„ (r) |') + r?|AB„ (r) |') dr < 2K,+^^hn-^ = K, + k2 is). 

(3.18) 
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3.3 i7^-Estimate of Um 

Wc establish a uniform upper bound for the iJ'^-norni of it„j by providing the estimate for the vorticity 
q„ = V X Vm- The Galerkin approximation (3.2a) is equivalent to 

+ vAVm - Pm {Um X Qm) = PmB {Bm, -Bm) • 

Taking the curl of the above equation we obtain 

^ + uAq^ -WxPm{umXqJ = Vx P^B {B^, B^) • (3.19) 
We use that in periodic boundary conditions 

/ {W X (j)) ■ ij)dx = I ^-{V X tp)dx (3.20) 

and for divergence free vectors 

V X (0 X V) = -(0 • V)V' + (V' • V)0. (3.21) 

Taking the inner product of (3.19) with q^, using that V • q^ = and the identities (3.20), (3.21) and 
(2.6), we reach 

ljt\lm\^ + "Wlmf = {B (g„, Um) , qj + (B {Bm, B^) ,V X qj . 
We bound the right hand side using (2.4), Young's inequality and (3.6) 

\{B {q^,Um) ,qm)\ < c\qJ'^^\M\\\qJ\^/^ 

<cu-^a-'kj\qj' + ^\\qj\' 

and by (2.9) 

|(B(B„,B„),V X qj\ < c||B„|p/2|AB„|V2||g^|| (3.22) 

<ct.-i|AB„|2 + ||B„f + ^||q„f. 

Note that since V • Vm = and due to the periodic boundary conditions we have 

\qm\ = |V X Vm\ = |VUm| = \\Vm\\, 

llGIlCG 

IqJ' < \\um + a^Aumf < (Ar^ + a^f \A'/^Um\'. (3.23) 

Hence we obtain 

~\qj' + '^hmf < cy-'a-'kl (Af^ + a'f \A'''um\' + cz.-^|^B„|2 + WB^f. (3.24) 

In the following we denote by Cj some constants depending on j/, 77, a, /ci, Ai. Integrating over {s,t) and 
using (3.18) and (3.17) we have 

\lm it) r < \Qm (^) r + Co (2K^ + ^fci/x-') + 2 (i^i + ^fciM"') dr. (3.25) 
We integrate this expression with respect to s over (|,i), f > and use (3.18), (3.23) to obtain 

|9„^WP<^^?i + Ci + y + | (3.26) 
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For t > ^ we integrate (3.25) with respect to s over the interval (t — -[;j^,tj. Note that, by applying 
also (3.18) and (3.23), we have 

it) l^<c,+c,(t-±-y+ celn (l - ^) ' • (3.27) 

From (3.26) and (3.27) we have, for t > 0, 

Wm it) I' < h it) , (3.28) 

where (t) has the foUowing properties 

(i) ks (t) is finite for all t > 0; 

(ii) ks (t) is independent of m; 

(iii) If either m^" ^ D (A^/'^) or B'" ^ V, then ks{t) depends on !/,r?,Q!, ||ix'"||, |B™| and 
limt^o+ ^3 it) = oo: 

(iv) limsupj^oQ ks {t) = R"^ < oo, E? depends on u, r], a, but not on u'" and B*". 
Returning to (3.24) and integrating over {t,t + t) , for t > 0, t > and using (3.28) we obtain 

y \\qj? <ki{t,T), (3.29) 
where kA{t,T) as a function of t satisfies properties (i)-(iii) as ks{t) above. 

Remark 3.2. If B*" G V and ti'" e £'(^), then by (3.16), Young's and Poincare inequalities we can 
bound (3.22) by 

\{B{B^,B^),V X qj\ < cu-'X-'/'k2\AB^\' + '^\\qj\'. 

Hence we have 

+ '^hmf < cv-'a-'kl {\-,'+a^f \A'"urn? + cv-^X-,"\,\AB^\' 
and by integrating over (0,t) and using (3.16) we obtain 

{im (t) I' < iQm (0) I' + cu-'a-'^kl (Ar' + a^)' (K^ + k^) + cv-^X^^'\2r^ [K^ + k^) . 
If, additionally, G D (A^/^), then using (3.23), we obtain 

Wm (t) I' < (Ar' + a')' + ci^-^a-''kj (Af' + q')' (ifi + fcz) + ci^-^X^^^^k2V^^ (i^i + fe) . 

(3.30) 

3.4 Existence of weak solutions 

Let us summarize our estimates. For any T > we have 

(i) From (3.6) 

ll^"i|lL°°([0,T];i/) - ^1' ll^™llL°°([0,T];y) - ^ °^ ll'"'"llL~([0,T];y') — ^"^^^ + a^) , (3.31) 

||-Bm||i,oo([o,T];H) ^ '^l- (3.32) 
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(ii) Prom (3.7) we have 



ll^mlll,2([0,T];y) - ^ 



\\u„ 



< 



or 



and 



IL2([0,T];D(A)) ^ 
fcl 



< 



1 I „2\ 



lL2([0,T];i/) ^ (Ai^ +Q;- 

ki 



< 



^\l^[0,T];V) - 



(3.33) 

(3.34) 

(3.35) 
(3.36) 



(iii) From (3.17) we have for any r e (0, T] 



Il~([t,T1;D(A)) 



< 



k2 ir) 



or Vr, 



< 



ka (r) 



and 



||S„|lioo([^_r];V) ^ ^2 (r), 



where ^2 (r) — > oo as r ^ 0"*". 



Now we estabhsh uniform estimates, in to, for Let us recall (3.2a). We have, by (3.35), 



Also, by (2.15), 



PrnB (Um^Vm) <c(Ai) || l^n 

D(A)' 



hence, applying (3.31) and (3.35), 



Pm-B (^mi ^m) 



L2([0,T];n(A)') 



< c 



1/2 



Additionally, by (2.13), we have 

<c{\i)-^'^\Bm\\\B^ 

therefore, using (3.32) and (3.36), we obtain 

2 kl 

\\PmB (-Bm,,-Bm)||L2([o^T];Z)(A)') ^ '^"7172 ' 



Consequently, by (3.2a) and the above 



dVr 



dt 



i2([0,T];Z)(A)') 



a^i/Ai 



1/2 



2a2 



7?Ai 



1/2 



and, in particular. 



dUri 



dt 



< 



L2([0,T];H) 
dB„ 



K 



(3.37) 



(3.38) 



Now we establish uniform estimates, in to, for . Let us recall (3.2b). We have, by (3.36), 

kl 



||^-Bm||j^2([o,T];y') ^ ^ 
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Also, by (2.4), 

\\PmB {Um, Bm)\\v' < C (Ai)"'/^ \\Um\\ ||-B„||, 

Hence, by (3.31) and (3.36), 

2 kf 

\\PmB {Um, Bm)\\L2no^T];V') ^ ^ 

Similarly 

\\PmB{Bm,Um)\\v' <c{Xi)-^/^\\Bm\\\\Um\ 

and 

2 kl 

\\PmB {Bm,Um)\\L2nQ^j'].Yi\ < C .1/2' 

Hence, from the above and (3.2b), we have 



dBr. 



dt 



2 u2 



From (3.34) and (3.38), using Aubin's Compactness Lemma (see, for example, [9, Lemma 8.4], [34] or 
[46] ) , we may assume that there exists a subsequence Um' of Um and u G ( [0, T] ; D [A) ) n C ( [0, T] ; ff ) 
such that 

Um' u weakly in ([0, T] ; D {A)) , (3.40a) 

Um' u strongly in ([0, T] ; V) and (3.40b) 

Um' u strongly in C ([0, T];H), (3.40c) 

as m' — ^ oo. Moreover, (d/dt) — > {d/dt)u weakly in {[0,T]; H). Or equivalently, by (3.35) and 
(3.37), there exists a subsequence Vm' of t>m such that 

Vm' V weakly in {[0,T]; H) , (3.41a) 

^ i; strongly in ([q, T] ; F') , (3.41b) 

Vm' ^ v strongly in C ([0, T] ; D (A)') , (3.41c) 

((i/dt) Vm' (d/dt) V weakly in ([0, T] ; D{Ay), as m' ^ oo, where v = u+a^Au is in ([q, T] ; i?)n 
C ([0, T] ; D {A)') . Also, by (3.36) and (3.39), there exists a subsequence Bm' of B„ and B e L'^ ([0, T] ; y)n 

C{[0,T];V') such that 

Bm' B weakly in ([q, T] ; F) , (3.42a) 
Bm' B strongly in ([0, T];H), (3.42b) 
Bm' B strongly in C ([0, T] ; F') (3.42c) 

and (d/dt) Bm' {d/dt) B weakly in ([q, T] ; F'). as to' ^ cx). 

Since — >■ 'U weakly in L^([0,T];iJ) and strongly in L^([0,T];y) and Bm' B weakly in 
([0, T] ; y) and strongly in ([0, T];H), then there exists a set ii^ C [0, T] of Lebesgue measure zero 

and a subsequence of Vm', Bm', which we relabel Vm, Bm respectively, such that Vm (s) v (s) weakly 

in H and strongly in V for every s G [0, T] \E, and Bm (s) ^ B (s) weakly in V and strongly in H for 

every s e [0, T] \^. 

Let w G D (A), ^ &V, then by taking the inner product of (3.2a) with w, and of (3.2b) with ^ and 
integrating over the interval [to,t] , t,to € [0,T], we have 

{v^ (i) , w) - {v^ {to) M + {b {u^ (s) , {s)) , Pmw) ds (3.43a) 

+ 1// {v^{s),Aw)ds= f {B{B^{s),B^{s)),Pmw)ds, 

J to J to 

(B„ {t) , ^) - {B^ (to) ,^)+ f {B (m„ (s) , B^ {s)) , Pmi) ds (3.43b) 

- f {B {B^ {s) , {s)) , P„^) ds + v f {{B^ is) ,^))ds = 0. 

J to J to 
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First we consider (3.43a). Since Vm (s) — ^ v (s) weakly in H, then for t, to G [0, T] \E 
(^m (0 ' ■"') ~ (^m (*o) ,w) ^ {v (t) ,w) — {v (to) ,w) , as m ^ 00 
and since w G D {A) we also have 



Now 



m— »oo 

For the nonlinear terms we have 



lim / (vm (s) , Aw) ds = (v (s) , Aw) ds. 



lim \PmAw — Aw\ = Km UPmif — = 1™ \PmW — w\ = 0. 



^ (S (m^ (s) , i;^ (s)) , Pmw) - (w (s) , V (s)) , ds 



Ho 

< 



+ 



B {u (s) ,v {s)),PmW — w) ds 

to > " " / D{Ay 

/ ( B (u(s) (s) — V (s)) ,w) ds 
Jto^ low 



=. /(I) + 7(2) . 7(3) 



By (2.15) 



4'^ <c(Ai)-'/' /V„ (s)IIK (s)||PmA«,-Aiz;|ds, 

"'io 



using Cauchy-Schwarz inequality we obtain 

I^^^ <c{\{)-^l^\P^Aw - Aw\\\u^ 



lL2([0,T];y) ll''mllL2([o,T];i?) 



hence by (3.33), (3.35) and (3.44) Hm, 
Again, by (2.15), 



■m — '*oo -*-m 



Ira — 'J* 



4'^ <c(Ai)-^/^ / \\u^{s)-u{s)\\\vrr.{s)\\Aw\ds 

J to 

and by Cauchy-Schwarz and (3.35), 



2m2 



T \ 

||it„(s) -m(s) f 



hence limm^oo Im) = 0, since ^ u in ([0, T] ; V). 



Finally, we show that lim. 



r(3) 



(3.44) 



■m— '■oo -^m 



4>{h) 



0. We define a hnear functional for h G {[0,T];H) by 



B{u(s),h),w) ds, 
to ~ /-D(^)' 



by (2.15) and Cauchy-Schwarz 

\<P {h)\<c (Ai)-i/* \Aw\\\u is) hmonv) \\h {s)\\ 
hence, due to (3.33), ^ is a bounded linear functional, and thus, since Vm — > v weakly in ([0,T] ;iJ), 

lim (j) {vm (s) -v{s)) = 0. 
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and hence limTO^oo Im) = 0. It remains to pass to the Umit in the right hand side element of (3.43a). 
/ {B {B^ (s) , (s)) , P^w) - {B {B (s) , B (s)) , w)y, ds 

J to 

< [ {B{B^{s),B^{s)),PmW-w)y,ds 

Jto 

[ {B{B^{s)-B{s),B^{s)),w)y,ds 



to 



{B{B{s),B^ {s)-B{s)),w)y, ds 
— 7(1) J. 7(2) I 7(3) 

Now, by (2.4) and Poincare inequaUty 

Ji'^ < c(Ai)-^/^ \\P^w - w\\ \\B^ ||i.([o,T];V) 

hence, by (3.44), lim„_oo Jm"* = 0. 
By (2.4) and Poincare inequahty 



Ji'^ < c(Ai)-^/^ /* \\B^ is) -B{s) \\\\B^ is) \\\\w\\ds 



and, applying Cauchy-Schwarz and (3.36), 



^S'<c(A,)-"ll»ll^(/ 



T \ 1/2 

\B^{s)-B{s)rds] . 



hence, since B^ — > B weakly in ([0, T] ; y), we have limm^oo Jm^ = (similarly to the argument 
given for 7^^). 



Similarly wc can show that limm^oc J-i 



(3) 



0. 



It remains to pass to the limit in (3.43b). Note that 

lim |P„$-^|=0. 



(3.45) 



We recall that B^ (s) — > B (s) weakly in V and strongly in H for every s G [0, T] \E, hence the 
convergence for the linear terms is easy. For the nonlinear terms we have 



(B (m„ (s) , B^ (s)) , P„^) - (B (u (s) , B (s)) , ^) ds 



to 
< 



[\b{u^ {s),B^ {s)),P^^-^)ds 

Jtn 



{B{u^ {s)-u{s),B^ {s)),Ods 



to 



+ 



[\b{u{s),B^ {s)-B{s)),0 ds 

Jto 



= . I ^(2) I 0(3) 

• ' m ' m ' 

Now, by (2.9) and Cauchy-Schwarz inequahty 

< c(Ai)-^/^ |P„^ - ^1 ||«„|L.([o,T];D(A)) 

hence, by (3.34), (3.36) and (3.45), lim^^oo = 0. 



\L^{[0,T];V) 
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Again, by (2.9) and Cauchy-Schwarz inequality 

5'^^ < c(Ai)-^/^|C| \\U^ (s) -«(s)lli2([o,T];D(A)) 

hence, since Um —>■ u weakly in ([0, T] ; D (A)), we have lim^^oo •S'm'' = (similarly to the case 
for I^). By similar arguments, using that B weakly in L^([0, T];!/), we obtain also that 

lim 'i'^^-' — n 

For the term J^^ {B {B^ (s) , (s)) , Pm^) we can perform the same estimates using the (2.10) to 
bound operator B. 

Hence, we can pass to the limit in (3.43) and we obtain that for every t, to & [0, T] \E 

{v (t) ,w)-{v (to) ,w)+ [ (b{u (s) , V (s)) ,w) ds + u [ {v (s) , Aw) ds (3.46a) 

Jto ^ ' -D(A)' Jt^ 

= f {B{B{s),B{s)),w)y,ds, 

Jto 

{B {t) ,^)-{B {to) ,0+ [ {B{u (s) , B (s)) ,Ods- [ {B {B (s) , u (s)) , ds (3.46b) 

Jto Jto 

+ T, f\{B{s),O)ds = 0. 

Jto 



for every w G D (A) ,^ gV. 

Now wc show that v eC ([0, T] ; V) (or equivalently m £ C ([0, T] ; V)) and B e C ([0, T];H). 

Notice that since ||'Um|lioo([o T]-y) — ^ {-^i^ + a^) and Vm v strongly in {[0,T] ;¥') then 
l|i'llz,~([o,T];V") ^ ^ (^r^ Hence (3.46a) implies that v (t) G Cy, {[0,T];V') because D {A) is 

dense in V. Since, also, for a fixed to, \\v {t)\\y, \\v (to)|lv"> as t — > to, then we have v € C ([0, T] ; V), 
or equivalently u G C {[0,T];V). 

Similarly, since ||-Bm||^oo(p,r];H) — ^^'^ — > B strongly in ([0,T] ; H) and D [A) is dense in 
H and because of (3.46b) we have B & C {[Q,T\-H). 

3.5 Uniqueness and continuous dependence of weeik; solutions on the initial 

data 

Next, we show the continuous dependence of weak solutions on the initial data and, in particular, the 
uniqueness of weak solutions. 

Let u, B and u, B be any two weak solutions of (2.17) on the interval [0,T] with initial val- 
ues u (0) = M^", B (0) = B™, u (0) = m"^ B (0) = B'". We denote v = u + a^Au, v = u + a^Au, 
6u = u — u, 5v — V — V and 5B — B — B. Then (2.17) implies 

j^5v + vA5v + B {5u, v) + B {u, Sv) = B {SB, B) + B {B, SB) , (3.47) 

^SB + r]ASB = -B {5u, B) - B {u, SB) + B {SB, u) + B {B, Su) , (3.48) 

(5ii(0) = ^w™ = (3.49) 

SB (0) = = - B". (3.50) 

Smcedv/dt€L^{[0,T];D{A)'), Su e {[0, T];D {A)) and dB/dt € L"^ {[0,T];V'), B, B, SB G 
L'^ ([0, T] ■ V) and due to the identities (2.6) and (2.8), we have for almost every t € [0, T] 

d 

dv, du ) +v 

I D(A)' 

B {SB, B) , Su)j,^^y + {B (B, SB) , Su)^, , 



Sv,Su) +v{\\Suf + a^\ASu\^) + (b{u,Sv),Su) 



j^SB, Sb"^^^ + r]\\SBf = -{B {Su, B) , SB)y, + {B {SB, u) , 5B)y, + {B (B, Su) , 5B)y, . 
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Notice that by theorem of interpolation by Lions and Magenes, see, e.g., [46, Chap. Ill, Lemma 1.2], 



■^5v, 5u . 

' D{A)' "'^ 



and 

thus we have 

d 



,\6u\'' + a^Wduf) + {\\6uf + a^\ASu\^) + (^B {u,Sv) ,Suj ^^^^^ 
= {B {SB, B) , 5u)^,j,y + {B {B, SB) , Su)^, , 



(3.51a) 



dt 



+ r]\\SB\\^ = -{B {Su, B) , SB)y, + {B {SB, u) , SB)y, + {B (B, Su) , SB) . (3.51b) 



By summation of (3.51a) and (3.51b) we obtain 
d 



dt 



{\Su\^ + a'^WSuf + \SB\^) + V {\Suf + c?\ASu\^) + ii\\SB\\ 



B{u,Sv),Su) +{B{SB,B),Su) + (B {SB,u) ,SB)y, . 

I D{A)' \ I D{A)' ^ 



By (2.16) we get 



B {u, Sv) , Su 



< c 



and by applying Young's inequality 



B {u, Sv) , Su 



< 



(\\u\\^'^\Au\^'^\Sv\\\Su\\ + \Au\\Sv\\Su\^'^\\Suf'^'^ , 



lAtlp (l^np + a-'WSuf) + -\\Sur + ja^l^-^nf 



1/2 



By (2.13) and Young's inequality we have 
' B{SB,B),Su\ 



<c\5Bf\\Bf + ^\\Suf + ^a'\ASu\. 



Also, by (2.10) and Young's inequality we obtain 



Summing up we have 



\{B{5B,u),SB)\ < ^\Au\'\SB\^ + ^\\SBr 



I {\Suf + a^Suf + \SBf) + ^ iWSuf + a'\ASuf) + ^\\SB\\ 



< 



-^\Au\^ + 1 + c||B||^ + -^\Au\A {\Suf + a'WSuf + \SBf) . 



We denote 



z{s) 



1/2 



1/2 



\Au\ 



and use Gronwall's inequality to obtain 

\Su{t)f +a^\\Su{t)f + \SB{t)f < {\Su{0)f +a^\\Su{0)f + \SB{0)f)exp(^J^ z{s)dsj , (3.52) 

since u,u G L'^ ([0,T] ; D {A)) and B € L'^ ([0,T] ; V) the integral J* z (s) ds^ is finite. Hence (3.52) 

implies the continuous dependence of the weak solutions of (2.17) on the initial data in any bounded 
interval of time [0,T]. In particular, the solutions are unique. 
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3.6 Strong solutions 



Theorem 3.3. Let T > 0, li™ e V, B'" e H. Then there exists a unique solution u,B of (2.17) on 

[0,T] satisfying 

u e i^, ((0, T] ; DiA^')) n LL ((0, T] ; n C ([0, T];V)n ([0, T] ; D{A)) (3.53) 

anrf 

B e L,- ((0, T] ; t/) n LL ((0, T] : n C ([0, T] ; il) n ([0, T] ; V) . (3.54) 
// B" s y and € -D(A) tften i/ie solution is the strong solution 

u€C ([0, T] ; i^(A)) n L\[0, T] ; /^(A^/^)), 
BeC([0,T];T/)nL2([o,T];i?(A)). 
//, additionally, u*" e D(A3/2^ ^f^^^ 

u e C([0, T] ; £)(A3/2)) n ([0, T] ; ^'(A^)) . 

Remark 3.4. Following the techniques presented in [16] (see also [13]) we can show that for any t > the 
solution is analytic in time with values in a Gevrey class of regularity of spatial analytic functions. As a 
result, we have an exponentially fast convergence in the wave number m, as m — > oo, in a certain sense, 
of the Galerkin approximation to the imique strong solution of (2.17), see, for instance, [10,27]. This 
Gevrey regularity result also implies the exponential decay of large wavenumber modes in the dissipation 
range of turbulent flows [11]. 

Proof. We use the Galerkin estimates derived in the previous subsections and similar ideas and compact- 
ness theorems in the corresponding spaces to converge to the strong solution. For (3.53) and (3.54) we 
need the estimates (3.28), (3.29), (3.6), (3.7) and (3.17), (3.18), (3.6), (3.7). For B'" € V, G D{A) 
we use the estimate (3.16) and if w*" G D{A^/^) we use (3.30). Also, since the strong solutions are weak, 
by uniqueness of weak solutions the strong solutions are unique. □ 



4 Convergence to the solutions of MHD equations as a ^ 0+ 

We emphasize again that our point of view is that the alpha model is to be considered as a regularizing 
numerical scheme. The next theorem shows that using the a priori estimates established previously, one 
can extract subsequences of the weak solutions of system (2.17), which converge, as a — > 0+, (in the 
appropriate sense defined in the theorem) to a Leray-Hopf weak solution of the three-dimensional MHD 
equations on any time interval [0,T]. For the definition and existence of weak solutions of the 3D MHD 
equations, see, for instance, [12] and [42] . The notion of a Leray-Hopf weak solution of MHD that satisfies 
the energy inequality (4.1) is inspired from a Leray-Hopf solution of NSE and formulated in the theorem. 
Also, if the initial data is smooth we prove that a subsequence of the strong solutions of the MHD-o; 
equations converges to the unique strong solution of the 3D MHD on an interval [O, T*(u™, B'")] , which 
is the interval of existence of the strong solution. 

Theorem 4.1. Let T > 0, ti*" G V, G H and denote by Ua, and — + a^Aua the 
weak solution of (2.17) on [0,T]. Then there are subsequences Uaj, Vaj, B^^ and a pair of functions 
v,B gL°° ([0, T];H)n ([0, T] ; V) such that, as Uj 0+, 

(i) Ua- — >■ V and Baj — > B weakly in L"^ ([0, T] ; V) and strongly in L"^ ([0, T];H), 

(a) — >■ V weakly in L^ ([0, T] ; H) and strongly in ([0, T] ; V) and 

(Hi) (t) V (t) and B^j (t) B (t) weakly in H and uniformly on [0,T]. 

Furthermore, the pair v,B is a Leray-Hopf weak solution of the MHD equations 

di3 ~ 

— + B{v,v)+iyAv = B{B,B), 
^ + B{v,B)-B {B, v) + r]AB = 



19 



with initial data v (0) = w'", B (0) = B'", which satisfies the energy inequality 

\v (t) |2 + \B (t) \'+2 f + v\\B{s)r) ds < \v (to) \' + \B (to) \' (4.1) 

J to 

for almost every to, <to <T and all t G [to,T]. 

Proof. From estimates (3.6) and (3.7), by passing to the limit (using the proof of Theorem 3.1), we have 
that the solution of (2.17) satisfies 



and 



{t) p + a^u^ {t) f + \B^ (i) p < fci 



2 r {i^{\M)f + a^\Au^{t)f) + v\\Bc,{t)f) dt < ku 
Jo 



notice that since a ^ 0+ we can assume that < a < L; consequently, we can bound the right hand side 
by fci := + + which is independent of a, therefore we can extract subsequences 



, , Baj , such that 



Uaj u weakly in {[0,T];V) , 
v„ . V weakly in ([0, T] ; H) and 
Ba, ^ B weakly in ([0, T]-V) , 



as Uj 0+. 



Now we establish uniform estimates, independent of a, for dB^/dt and dua/dt. From (2.17b) we 
have 

_idB^ 



notice that by (2.13) 
hence 



and similarly 

Hence 



< \\A-^B{u„,B^)\\ + \\A-^B{B^,u^)\\+rj\\Bj, 
|A-iB(u„,B„)| <cAr^/Vall|Sa||, 

|-B(Ma,-Ba)||i2([o,T];Z)(A)') <C>'1^'^ \Uo,{t)\^\\Bo,{t)\\'^dt 

dBr 



dt 



<K, 



L^{[0,T];D(A)') 



where K is independent of a. 
From (2.17a) we have 



j^-ldUa 



dt 



< 



A-^{l + a^A) ^B{u,v] 



+ 1^ uJ + 



A-^{l + a'^A) ^B{B,B] 



and using (2.15) 



\A-^I + a^A) 'B{uc,Va)\< \A-^B{uo,,v^)\ 

< c\'^^^^\\Ua\\\Va\ 

< cX^^^^\\Ua\\ {\Ua\ + a^lAUal) , 



thus 



\A-' (/ + a' A) ' B (««, v^) |2 < 2cA-^/'fci {\\u^f + a^\Au^\'') , 
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and 



J 

Jo 



\A-^ {I + a'^A) B {uc {t) , (i)) < cX^ ' k^iy~^. 

JO 

Also by (2.13) 

As a result we have 

du^ < K. 

LHIO.TYMA)') 



dt 



Using Aubin's Compactness Lemma (see, for example, [9, Lemma 8.4]) we can extract subsequences of 
and , which we relabel by Ua- and Baj respectively, such that — >■ u and Bcj — > B strongly 
in L2 ([0, T] ■ H) and strongly in C ([0, T] ; D {A)') , as aj 0+. 
Observing that 

\\Va, - Ma, \\l2{[0,T];V') = H^"^ W H'^* - "j^^^^^l' 

we obtain that v^^ ^ w in ([0, T] ; V), as aj — > 0+; and hence also that u{t) = v{t) almost 
everywhere on [0,r]. 

Now, following the lines of the proof of Theorem 3.1, we can extract further subsequences (which we 
relabel by Uaj,Vaj and B^j) and show that as aj — > 0+, 

B {Uaj ,Va^) ^ B{V,V) = B {v, v) 

weakly in ([0, T] ; D (A)') , and 

B{Bc,.,Ba,) ^ B{B,B), B{ua,,Bc,.) ^ B {v,B) , B {Bc,.,Uaj) -^B{B,v) 
weakly in ([0,T] ;V'). Hence, we can pass to the limit (in the interpretation given by (2.19)) in 

jVa, , '^)^^^^, + ' ) ' '^)^(^), + " (^«. ' = (-^"^ ' ) ' '^)y' ' 

, + {B {ua, , B„J , ^) - (B (B„ . , n„,) , ^) + r/ ((B„, , ^)) = 0, 
w G D (A) , ^ e V and we obtain that 

^v, w\ +{B {v, v) , w) jjr^y + V {{v, w)) = {B {B, B) , w)y, , 

/ D(A)' 

jB. i^^^ + {B {v, B),0- {B [B, v),0+rj {{B, 0) = 0, 

for every w £ D {A) , £V and for almost every t S [0, T]. 

Now, since v e Lp- {[Q,T]\V), one can show that B{v,v) e L^([0,T];F') and then also that 
{d/dt)v G {[0,T] ;V'), and since w G D{A), which is dense in V, we obtain the weak formulation 
of the MHD equations 

^v, + {B {v, v),w)y,+u {{v, w)) = {B [B, B) , w)y, , 

(j^B, i^^^ + {B {v, B) , - {B {B, v),i) + ri {{B, 0) = 0, 

for every w, ^ &V and for almost every t £ [0, T]. 

We notice, that every weak solution of (2.17) satisfies the energy equality (3.1) and hence the energy 
inequality (4.1) follows by passing to the liminf as a ^ 0+, using the fact that if Xa — > a; weakly in a 
Hilbert space X, then ||a;|| < liminf □ 
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Theorem 4.2. Let T > 0, e D{A), G V arirf denote 6j/ and t»a + a^^^j^ f/ie 

strong solution of (2.17) on [0,T]. Tften there exist T^, — T^{fl,iy,r],u^'"-,B^"), < T* < T, subsequences 
, i'qj , and a pair of functions v,B G L°° ([0, T*] ; y) n ([0, T*] ; -D(j4)) smc/i that, as aj — > 0+, 

(i) Ua - V and B^j B weakly in ([0, T*] ; D{A)) and strongly in ([0, T*] ; V), 

(ii) Vaj V weakly in ([0, T*] ; V) and strongly in ([0, T*] ; H) and 

(Hi) Ua - (t) —>■ V (t) and Bq,. (t) B (t) weakly in V and uniformly on [0,T,]. 

Furthermore, the pair v,B is the unique strong solution of the 3D MHD equations on [0,T,] with 
initial data v (0) = it™, B (0) = S*". The strong solution of the 3D MHD equations satisfies the energy 
equality 

\v {t) p + \B (t) P + 2 /* {u\\v{s)f + v\\B{s)f) ds = \v {to) P + \B {to) |^ 0<to<t< T,. 

J to 

Proof. To prove the theorem we need to show that there exists T* such that we have a uniform (inde- 
pendent of a) bound on 

\\ua{t)f+a'\Aua{t)\^ + \\Ba{t)f (4.2) 

and 

(v{\Aua{t)\'' + a2|A3/2-ix„(f)|2) + n\ABa{t)\'''^ dt (4.3) 

in [0, T*]. Then we can continue similarly to the proof of the previous theorem, appropriately smoothing 
the data and replacing T by T*. Next we derive the formal estimates on (4.2) and (4.3) that can be 
proved rigorously using the Galerkin approximation scheme and then passing to the limit using the proof 
of Theorem 3.1. 

Let us recall (3.13). By (2.9) and several applications of Young's inequality we bound 

(i3(tx«,fa),^Ma) < c\\ua\\^''^\Aua\^/'^{\\ua\\ + a^\A^''^Ua\)\Aua\ 

< cv-^\\uat + v-''a'\Au^\^ + + '^a^\A^/^u^\\ 

By (2.9) 

\{B {B^,B^),Auo:)\<c\\BX'^\AB„\^/^\\Ba\\\Aua\ 
By (2.9) we also have 

|(B (B,, tx„) , AB„)| < c||B„ f /2 II I |3/2 

<cr?-3||B,||6+j?-3||tx„f + ||AB„|2 

and by (2.10) 



|(B(m„,B„),AB„)| <c||B„f/2||tx,|||^B,|3/2. 
Hence by (3.13) and the above estimates we have 

^ (llixaf + a^l^Wal" + ||B„||2)+z. {\Auc,\^ + a2|A3/V|')+r?|AB„|' < c/x"^ (||u„f + a^\Au^f + \\Bo 

(4.4) 

Denote 

y=\\ua\\^ + a^\Aua\'' + \\Baf. 

Now, if y(0) = 0, that is m'" = B*" = 0, then the solution is steady Ua{t) = 0, Va{t) = 0, Ba{t) = 
and v{t) = 0, B{t) = exists for all t>0. Otherwise, from (4.4) we have 

d _o o 



22 



and thus 

y{t)<2y{0) 

for < t < fcV^y (0)"^ We conclude that 

\\ua (t) f + a^\Auc, + IIBaWf < 2 iWu'^f + a^l^u^l^ + ||B-f ) 
for < t < T* := min (t, |c^/i^j/ (0)~^^ . Also, by integrating (4.4) over (0,T*), we obtain 

1^' (u{\Aua{t)f + a^\A'/^u„{t)f) + v\AB^{t)f) dt 
Assuming that < a < L, the bounds are independent of a. □ 

5 Discussion 

We proved the well-posedness of the three-dimensional MHD-a model (1.6) in the periodic boundary 
conditions. This model modifies the nonlinearity of the MHD equations (1.1) without enhancing dissipa- 
tion. Wc showed that the model has a unique global weak (or strong, for smooth initial data) solution. 
Also, there is a subsequence of weak solutions of the MHD-a equations that converge, as a ^ 0+, (in 
the appropriate sense) to a Leray-Hopf weak solution (which satisfies the energy inequality (4.1)) of the 
MHD equations (1.1) on any time interval [0,T]. Also, if the initial data is smooth, a subsequence of 
solutions converges for a short interval of time, to the unique strong solution of the MHD equations on 
this interval. These properties are essential for the a models to be regarded as regularizing numerical 
schemes. In a follow up paper, wc intend to do the error estimates in which we will investigate the error 
in terms of to and a. Namely, the distance between the solution of the Galerkin MHD-a model to that 
of the exact strong solution of the MHD equations, for smooth initial data. 

There arc many different a models. For example, the global wcll-poscdness can be shown for the 3D 
Modified-Leray-a-MHD model (1.9). However, at the moment we are unable to find a conserved quantity 
in the ideal version of (1.9), which can be identified with a cross helicity, contrary to the MHD-a model 
(1.6), where there exist the ideal invariants that could be identified with the three invariants of the 
original 3D MHD equations. 
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